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We formulate a solvable model that describes generation and saturation of mean magnetic field in 
a dynamo with kinetic helicity, in the limit of large magnetic Prandtl number. This model is based 
on the assumption that the stochastic part of the velocity field is Gaussian and white in time (the 
Kazantsev-Kraichnan ensemble) , while the regular part describing the back reaction of the magnetic 
field is chosen from balancing the viscous and Lorentz stresses in the MHD Navier- Stokes equation. 
The model provides an analytical explanation for previously obtained numerical results. 



o 
o 



(N 



> 

m 
in 
o 



^ ■ 

Or 

6 ■ 
a: 



1. Turbulent motion of a conducting fluid generates 
magnetic fields in planets, stars, and interstellar medium. 
Observational |jj and theoretical data suggest that 
simultaneously with developing intense small-scale fluc- 
tuations, magnetic fields appear to be ordered on scales 
much larger than the correlation length of the turbulent 
velocity field. Explanation of the large-scale field has 
remained a challenge for the astrophysical dynamo the- 
ory for a long time. A leading approach suggests that 
the large-scale field is generated due to broken parity 
in the turbulent velocity fluctuations. Mathematically, 
this means that the fluctuating velocity field u possesses 
nonzero helicity, i.e., J u-(Vxit) ^ 0. The analytical con- 
firmation of the large-scale field generation is based on 
the splitting of the magnetic field into the large scale and 
small scale components (compared to the velocity corre- 
lation length) , and averaging the induction equation over 
the small scales || . This procedure leads to the equation 
for the mean magnetic field {B l ), which can be written 
in the general form: 



3 t (B')=Vx (a(B l ))+l3A(B t ) + .. 



(1) 



Functions a and (3 should be obtained from the theory. 
As an estimate, one can accept a ~ u and (3 ~ ul, where 
I is the characteristic scale of the velocity field. It is as- 
sumed that the mean field is changing slowly in space, 
therefore the higher order space derivatives of (B l ) are 
omitted in ([!]). The crucial fact is that the coefficient a 
describing the amplification of the mean field is zero if 
the velocity fluctuations possess no helicity M . 

There is no general theory that would allow one to 
derive coefficients a and (3 from the microscopic equa- 
tions. A self-consistent derivation of the mean field equa- 
tion can be carried out under certain simplifying assump- 
tions. One of them is the assumption about the geome- 
try of the problem (velocity field is assumed to be two- 
dimensional) , which leads to a solvable model introduced 
by Vainshtein || . Another assumption is about smallness 
of magnetic fluctuations compared to the mean field (see, 



e.g., ]10| , |1 lfl ). However, analytical and numerical results 
on the initial (kinematic) stage of the dynamo demon- 
strate that it is the small scale magnetic fluctuations that 
are most energetic, not the large-scale mean field ^H. 
Therefore, an analytical treatment of this opposite limit 
in a general 3-dimensional case is in demand. 

Moreover, the previous analytical and numerical inves- 
tigations raised the question of the so-called catastrophic 
a-quenching in the mean-field dynamo. The mean mag- 
netic field was shown to saturate when its energy reaches 
the equipartition value divided by the magnetic Reynolds 
number, R m For the critical discussion of the a- 

quenching see also [Q. The magnetic resistivity thus 
played an essential role in the saturation mechanism. As- 
trophysical applications, e.g., galactic dynamos, provide 
a unique physical setting of extremely large magnetic 
Prandtl numbers (ratio of fluid viscosity to magnetic dif- 
fusivity), reaching 10 14 -10 22 , where the back reaction of 
the growing magnetic field can come into play before the 
magnetic field reaches the diffusive scales. The defini- 
tive numerical investigation of this problem requires very 
high resolution of the energy containing fine scales and is 
presently without the reach. An analytical understand- 
ing of this limit is therefore important. 

In the present paper we introduce a solvable model de- 
scribing the nonlinear mean field dynamo in the limit of 
a large magnetic Prandtl number. We do not resort to 
any artificial scale separation procedure or to the special 
geometry, and do not assume that the magnetic fluctu- 
ations are small. Instead, we make simplifying assump- 
tions about the velocity field. We assume that the micro- 
scopic fluctuations of the velocity field are homogeneous, 
Gaussian, and short-time correlated, which constitutes 
the so-called Kazantsev-Kraichnan ensemble while 
the the regular, "back reaction," part is chosen from bal- 
ancing the viscous stress and the magnetic field stress. 
We also assume that the random velocity field possesses 
helicity. We believe that this model can provide some in- 
sight into understanding of the mean-field dynamo, and 
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can be complimentary to the previously developed ana- 
lytical treatments. 

We show that in our model the fluctuations of the 
magnitude of the magnetic field are independent of the 
fluctuations of its direction. In the presence of helicity, 
the probability distribution function (PDF) of the di- 
rection of the magnetic field becomes nonstationary and 
anisotropic in x space. However, the probability distri- 
bution function of the magnitude of the magnetic field 
does not depend on x. Although the anisotropic part 
of the PDF decays in the course of time, the magnitude 
of the magnetic field increases. The interplay of these 
two factors leads to the generation of the non-decaying 
mean magnetic field, (B l ). This exact mechanism of the 
mean magnetic field generation is the first main result 
of the paper. The mean magnetic field is shown to be 
determined by the geometry of the problem and by the 
boundary conditions; the fact previously emphasized by 
Blackman and Field 

When the Lorentz force becomes strong enough, one 
cannot neglect its back reaction on the velocity field. If 
the velocity field helicity is small, this back reaction af- 
fects only the probability density function of the mag- 
netic field magnitude, not the PDF of the magnetic field 
direction. As an example, we present a model based on 
the balance of viscous stress and magnetic stress, that 
can be solved exactly and clearly demonstrates this sat- 
uration: anisotropy of the PDF decays as before but the 
second factor, the magnetic field magnitude, saturates 
and cannot support the mean magnetic field anymore. 
The a term thus becomes reduced. This mechanism of 
saturation of the growing mean magnetic field is the sec- 
ond main result of the paper. 

In Section 2 we formulate the model and proceed with 
the detailed derivation of the mean field equation. In 
Section 3 we suggest a simple closure describing dynamo 
saturation within the present model. Section 4 discusses 
possible generalizations of the model and its limitations. 

2. Let us first neglect the Lorentz force and assume 
that the fluctuating velocity field u z (x,t) is random, ho- 
mogeneous, Gaussian, and short-time correlated: 

(u l {x, t)u k (x't')) = S(t - t')n lk (x - x'), (2) 

where due to the lack of mirror invariance we have K lk (x— 
x') 7^ n lk (x' — x). Magnetic field obeys the induction 
equation: 

B k u\ 



d t B l 



u k Bl 



nAB\ 



(3) 



where lower indices denote derivatives with respect to the 
corresponding spatial coordinates, and we assume sum- 
mation over the repeated indices. In the case of large 
magnetic Prandtl numbers (which is usually the case in 
astrophysical settings and is considered here), we neglect 
the magnetic diffusivity n. Using equation (|J) as the 
Langevin equation for magnetic field fluctuations, one 
can easily derive the Fokker-Plank equation for the prob- 
ability distribution function of the magnetic field mea- 
sured at some point x, P(B l ; x, t): 



d t P = n AP + k 2 B 2 



d 2 P 



d{B 1 ) 2 d + 1 



i nfe 



K 2 B B 



d 



2ge B % -—^SIiP. 



d 2 P 
dB l dB k 

(4) 



The coefficients in this equation can be read off from 
the following expansion of the velocity correlation func- 
tion (|^) in small argument y = x — x' for incompressible 
velocity field: 



«*(y) = K0 S ik - | (y 2 6* k - 2 -{^ 



+ ge M y l + ..., (5) 



where d is the space dimension (d — 3), and the g term 
describes kinetic helicity. Note that the value of g should 
obey the realizability condition 



9 < o K 0«2, 



(6) 



that follows from the fact that ((ou + V x u) 2 } > for 
any a. 

For completeness, we would like to present here the 
main steps of the derivation of Eq. (||) from Eq. (|^). 
First, let us introduce the so-called Z function defined 
as Z(X,x,t) = exp(iX l B z (x,t)). Obviously, the aver- 
age of this function over the random velocity field, it, 
or, equivalcntly, over the resulting magnetic field distri- 
bution P(B;x,t), gives the Fourier transform of the dis- 
tribution P(B;x,i) with respect to B l . Due to Eq. 
the Z function satisfied the following differential equa- 
tion: 



d t Z = -iXWBiZ + iX l u\.B k Z 



d 



d 



- u ^ z + uiy W Z - 



(J) 



Let us now average this equation with respect to the 
random field u. To do this we first formally iterate the 
Eq. (R) once, getting: 



d t Z = u k {t) 



dx k 











tW^(*0-«KOAW(O 



dx r 



dX k 



d 



, (8) 



where the integration is performed with respect to t' and 
we do not write the dependence on x in the right hand 
side of (||) . The average over u can now be done indepen- 
dently of the Z function, since the Z function depends on 
the velocities taken at earlier times, t' < t, and due to 
causality, cannot depend on u(t). We thus average u's 
and Z's independently in (g), and Fourier-transforming 
the resulting equation with respect to A, arrive at Eq. (Q). 
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Now, introducing new variables, the magnitude of the 



magnetic field B and its direction n l 
we cast the Eq. (^) into the form: 



B l /B (n 2 



d t P = k AP + L B P + L n P, 



where 



t 1 d d+1 d 

hB - K2 d+lB^dB B 8B 

a 2 



i) 



(9) 



(10) 



on l anr an 1 



(11) 



When differentiating with respect to n z , we assume that 
all components of n l are independent. This is a legiti- 
mate procedure since one can look for the solution in the 
form P = PS (I — n 2 ), and observe that the 5 function 
factors out, i.e., L n S(l — n 2 )f = 5(1 — n 2 )L n f for an 
arbitrary function /. 

The statistics of magnetic field magnitude and mag- 
netic field direction are independent if they are initially 
independent. Let us thus look for the solution of the 
Fokker-Planck equation (|j) in the factorized form 

P = P B {B,t)G(n,x,t), (12) 

where P B is the function satisfying the equation 

d t P B = L B P B . (13) 

This equation can be solved exactly to give the well- 
known log-normal distribution of B (see e.g. Jl^]), but 
we do not need this solution for our present purposes. 
Let us concentrate on the equation for the G function 



d t G = K AG + L n G. 



(14) 



If the kinetic helicity is nonzero, the last term in the 
operator ( [Tl| ) leads to coupling of n and x. A rigorous 
analysis would require solving the eigenvalue problem for 
the operator in the right hand side of (^|), which we are 
going to do elsewhere. For now, we can accept that the 
G function in ( |l2"| ) corresponds to the largest eigenvalue 
of the operator in the right hand side of (Q). The rele- 
vant structure of the solution can be simply understood 
in the following way. Let us first assume that G is inde- 
pendent of x. One can then formally expand G in powers 
of n l : 



G(n, t) = l + SV + Q lk n l n k + D lkl n l n k n l + 



(15) 



where S, Q, and D are functions of t. Plugging this 
expansion in the Fokker-Planck equation (|l4|), we easily 
check that the higher the order of the expansion coeffi- 
cient, the faster it decays, and therefore for large times 
we can safely truncate the expansion ( |l5| ) after the first 
relevant term, G = 1 + S' l n l . Higher order terms can how- 
ever be important for consideration of higher moments of 



the magnetic field but are irrelevant for our consideration 
of the mean magnetic field (B 1 ). This truncated solution, 
describing anisotropy of the magnetic field distribution, 
behaves as 



S oc exp[— K 2 (ei — 1)^]- 



(16) 



This suggests to look for the general solution of (|TJ) in 
the form 

G(n,x,t) = 1 + B l {x,t)n t exp[-K 2 (d - 1)*]. (17) 

Equation (|l|) now reduces to (in three dimensions) 

d t B l = k AB 1 + 2. 9 V x B\ (18) 

and we recover the equation for the mean field (|l]), with 
the eddy diffusivity j3 = kq, and a — 2g. 

To reveal the physical sense of the function B % (x, t), let 
us calculate the mean magnetic field (B l ) using the distri- 
bution (|l2|). At first sight, it should decay since the mag- 
netic field distribution function becomes more and more 
isotropic in the course of time, and the isotropic distribu- 
tion has no preferred direction for the magnetic field B l . 
However, the amplitude of the magnetic field grows mak- 
ing such an average nonvanishing. Indeed, multiplying 
the equation (|l3| ) by B d (magnitude B times the volume 
factor B d_1 ) and integrating with respect to B we get 
the average of the magnitude 



(B) = B exp[K2(d-l)t], 



(19) 



where we assumed (B) — B at the initial moment. Fi- 
nally, averaging with respect to angles with the aid of the 
distribution G, we obtain (B l ) = B l (x,t) B /d. Note 
that the exponential factors cancel out exactly in ( |l6| ) 
and (|l|). 

To conclude this section, we would like to explain the 
meaning of the averaging procedure in the above formu- 
lae. Mathematically, we average over the Gaussian en- 
semble (|2j). Physically, the averaging is understood to 
be performed over the space. More precisely, we average 
over the scales much larger than the velocity correlation 
length, but much smaller than the scales at which the 
function B l (x,t) is changing. Since on these scales the 
velocity field is not correlated, we are effectively averag- 
ing over non-correlated regions of space, which is equiv- 
alent to the ensemble averaging. This is a sensible phys- 
ical procedure, since, for example, the scale of the mean 
galactic magnetic field is > 1 kpc, while the velocities are 
correlated at the scales < 100 pc. 

3. Equation (|l8|) shows that the growth of the mean 
magnetic field is determined by the eigenvalues of the 
operator on the right hand side of (|l8|). The solution 
depends on the magnitude of the helicity of the velocity 
field, on the geometry of the problem, and on the bound- 
ary conditions. Space Fourier transform of the Eq. dlg| ) 
shows that the eigenvalues of the r.h.s. operator are: 
Ai = — K fc 2 ,A2,3 = —nok 2 ± 2gk. Choosing g to be 
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positive, we obtain that the largest growth rate is equal 
to Amax = <? 2 / K o an d is achieved at k = g/Ko if this k 
is allowed in the system. Obviously, this growth can 
not proceed for infinitely long time, since eventually the 
growing magnetic field will be brought to equipartition 
with the velocity field. We therefore need to introduce 
the back reaction of the magnetic field to the velocity 
field. We do not know any rigorous way of including this 
reaction into the Eq. (Q). A simple physical model can 
be suggested in the case when the mean magnetic field is 
much smaller than the magnetic fluctuations. This can 
happen when the following two conditions are met: 

(1) The helicity of the velocity field is small, 
i.e., kq/1,K21 3> g, where I is the correlation length of 
the velocity field. In this case, the generated mean mag- 
netic field is concentrated on the scales much larger than 
the velocity correlation length, and is growing with the 
rate much smaller than that of the fluctuations growth. 

(2) The initial mean magnetic field is smaller than 
the initial magnetic fluctuations. The realizability con- 
dition (|]) then ensures that the mean magnetic field is 
much smaller than the magnetic fluctuations at all con- 
secutive times. 

We now show that under these assumptions the 
equipartition is first achieved locally between the mag- 
netic fluctuations and the velocity field. In a Lagrangian 
frame, magnetic field is rotated and stretched by the fluc- 
tuating velocity gradient matrix, u\: 



dynamical stress 



dt 



u\.B k 



(20) 



To model the back reaction, we want to represent this 
fluctuating gradient as some random part, u\ plus some 
regular part (that does not include the random vari- 
able u\.) which is proportional to the magnetic fluctu- 
ations B l . Field B l , being much smaller than the fluc- 
tuations, would give higher order terms in g/^l, gl/^o, 
and should not be relevant at least at the beginning of the 
saturation. We have the following traceless tensors at our 
disposal: e ikl B l , B\, and T lk = B l B k - S lk B 2 /d, which 
are constructed from the local magnetic field. The first 
two tensors can not be physically relevant, since they do 
not contribute to the energy exchange between the mag- 
netic and velocity fields, which due to ([20j) is proportional 
to (B l B k u\). As we show below, the choice of T %k leaves 
the operator L n intact and changes only the Lb part of 
the Fokker-Planck equation (Q), (g). But before we pro- 
ceed with the derivation we would like to give a more 
sound physical motivation for choosing the regular part 
of u\ in the form of T tk , as suggested in [ p"7[ . 

Let us assume that the dynamo saturation starts when 
the viscous stress vAu 1 in the right hand side of the in- 
compressible MHD Navier-Stokes equation balances the 



V fe \B*B k --5 lk B 2 j -Vip. 

This balance seems to be reasonable at the onset of the 
back reaction, since it is the smallest fluid eddies that are 
less energetic and that are affected first by the growing 
magnetic field. Also, assuming that the pressure ensures 
the incompressibility of the velocity field (itl must be 
traceless) we are left with 



1 



B % B 



1 



8 lk B' 



'k ' 



(21) 



where the random part u\. obeys (||) , and the regular part 
is just ~ T lk . The structure of the regular part of the ve- 
locity gradient shows that the magnetic field is stretched 
along itself, but not rotated. It is instructive to see the 
possible physical meaning of this ansatz. Numerical re- 
sults of Jl8|-p0[ suggest that the growing magnetic field 
becomes organized in filaments that occupy small frac- 
tion of space. The length of the filaments is of the order 
of the size of energy containing eddies. The magnetic 
field amplitude inside the filaments is large, and the sat- 
uration with the velocity fluctuations occurs first inside 
these filaments. With such saturation, small-scale ve- 
locity fluctuations cannot bend the filaments, since this 
would lead to a large restoring Lorentz force. The reduc- 
tion of curvature of the filaments during the saturation 
process has indeed been observed numerically by Bran- 
denburg et al pl|. We thus have to again exclude the 
local rotation (terms e tkl B l and Bl) and leave only the 
stretching along the magnetic field lines which is given 

by©. 

As we have mentioned, the ansatz (^lj) does not change 
the L n part of the Fokker-Planck equation. The Lb part 
changes and takes the form: 



d t P B = L B P B 



d- 1 1 d 



B d+2 P B . 



(22) 



vd B d ~ x dB 

The saturated PDF can be found analytically as a sta- 
tionary solution of this equation, and turns out to be 
Gaussian, 



P B (B) = Acxp - 



-B' 



(23) 



where A is the normalization constant and d — 3, as 
oppose to the log-normal distribution at the initial, kine- 
matic, stage. Comparison with direct numerical simula- 
tions [Q shows that this gives the qualitatively correct 
behavior of the magnetic field PDF - the PDF indeed 
saturates at the Gaussian. 1 It is curious that the clo- 
sure (|2l|), suggested here for the onset of the back re- 
action, also works for the fully saturated regime. The 



1 A. Schekochihin was able to find the exact time- dependent 
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Gaussian stationary distribution of the magnetic field 
was earlier obtained in numerical simulations by Catta- 
neo pp[ | for the central part of the distribution function, 
while the tails of the PDF seemed to be non-Gaussian. 



We thus expect our local ansatz (21) to be valid only 
for the central part of the PDF; indeed, not being able 
to capture the spatial structure of the magnetic field it 
cannot account for the intermittency effects described by 
the non-Gaussian tails of the PDF. Note, that Gaussian- 
ity of the magnetic field distribution was also assumed 
in an analytical model for the mean magnetic field by 
Subramanian pl[ , see also numerics in p2fl . 

The mechanism of the saturation of helical dynamo 
with small mean field is thus the following. The growing 
fluctuations of the magnetic field magnitude, described 
by Ps{B,t), saturate locally with the velocity fluctua- 
tions, i.e., as described by the above model (|2l])-(22). Af- 
ter that the isotropisation of the magnetic field distribu- 
tion proceeds faster than the amplitude growth, and sup- 
presses further growth of (B l ): exponential decay of jiq ) 
is not compensated by (|19| ) anymore, and due to the re- 
alizability condition (p|), can not be compensated by the 
growth rate Amax either. This provides a mechanism for 
the observed reduction of the a factor. Remarkably, our 
analytical result based on ansatz (21) is also in accord 
with the physical pic ture of dynamo saturation advo- 
cated by Cattaneo |19| , po| : saturation occurs first locally 
in filaments that gradually fill the whole space. 

Finally, it is instructive to see at what level the mean 
magnetic field should saturate. Due to (|2^), the sat- 
urated value of magnetic fluctuations energy is esti- 
mated as the energy value of the smallest turbulent ed- 
dies, W ~ vk.2 ~ E/Re 1 ^ 2 , where E is the total fluid 
energy and Re is the Reynolds number of the fluid turbu- 
lence. The saturation time can then be found from ( |l9| ) 
to be t s ~ log(WyWo)/(4K 2 ). Assuming that the mean 
magnetic field grows with the maximal possible growth 
rate, Amax = 9 2 /kq, we obtain for the saturated energy 
of the mean field, W: 



W ~ Wo (W/W ) ' 



(24) 



where 7 = g 2 /(2k K2). Setting Wo ~ Wo, and letting 7 
have the maximal possible value allowed by the realiz- 
ability condition (|J), 7 m = 5/16, we estimate: 



W < E (Wo/S) 1-7 " 1 Re 



-7-/2 



(25) 



We would like to stress the strong dependence of the sat- 
urated mean magnetic field on the value of the initial 
magnetic field, and rather weak dependence on the fluid 
Reynolds number. Since in our model the back reaction 



starts before the resistive scales are reached, the answer 
is independent of the magnetic diffusivity. 

4. In the present paper we developed a frame- 
work for analyzing the generation and saturation of the 
mean magnetic field, which is based on the Kazantsev- 
Kraichnan ensemble for the velocity field. Supplemented 
by the physically motivated ansatz for the magnetic field 
back reaction, such a model accounts for the growing 
mean magnetic field exactly. 

Let us now discuss possible limitations of the approach. 
The solution of the mean-field equation requires the spec- 
ification of the particular geometry and can be different 
for different boundary conditions. Of course, the con- 
dition of homogeneity of the velocity fluctuations could 
be incompatible with real boundary conditions, which 
would force us to consider more general velocity correla- 
tors, k(x + x'; x — x 1 ) and to solve the eigenvalue problem 
analogous to (jl4|). Another generalization may be nec- 
essary for the case of arbitrary, not small, magnitude of 
the mean magnetic field. In this case, the ansatz for 
the regular part of the velocity field can include large B i 
field as well. Also, the neglected tensor B l kl though not 
leading to the energy transfer, can be important for the 
isotropisation of the PDF in this case. Another possible 
generalization is to consider a compressible velocity field. 
All these questions are for the future. 

I am indebted to Eric Blackman, Axel Brandenburg, 
Vladimir Paricv, and Dmitri Uzdcnsky for important 
suggestions and comments on both the physics and the 
style of the paper, and to Lars Bildsten and Alexander 
Schekochihin for the discussion of the results. 
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